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1. Introduction and summary of results 

Static non-linear susceptibilities serve to characterise a phase transition and its 
universality class, especially in spin-glass systems in which the linear counterpart does 
not show a divergence but just a cusp at a finite critical temperature. Experimental 
results for standard spin-glass samples have now been available for over 20 years pQ E] 
while studies of more exotic systems exhibiting spin-glass ordering, such as manganites, 
are currently being carried out [3]. 

Recently it has been pointed out that dynamic non-linear susceptibilities can also 
be very useful to characterise, and eventually understand, the slow dynamics of super- 
cooled liquids, their arrest [H EJ [7], and the non-equilibrium dynamics of the low 
temperature regime [SJ [9] . 

In order to separate the truly non-trivial behaviour of glassy systems (with 
or without quenched disorder) from phenomena also present in simpler cases, such 
as ferromagnetic domain growth and other simple phase ordering mechanisms, it 
is important to understand the behaviour of the static and dynamic generalised 
susceptibilities in solvable toy models. The aim of this paper is to present a detailed 
analysis of static and dynamic linear and non-linear magnetic susceptibilities in generic 
spherical disordered models with two-body interactions [10]- [SJ. We shall investigate how 
these quantities depend on the density of states of the two-body interaction matrix. In 
the study of susceptibilities, i.e. the influence of an external field h on the system, 
two thermodynamic limits in h and iV (the number of spins or system size) can be 
considered, 

(i) The applied field goes to zero first and the thermodynamic limit is taken after: 

lim lim . 

This is the limit which always coincides with static fluctuation dissipation relations 
relating susceptibilities to correlation functions. 

(ii) The field goes to zero after the thermodynamic limit is taken: 

lim lim . 

This is the limit of the critical isotherm. 

It will turn out that these two limits are equivalent in the high temperature phase, 
however they differ in the low temperature regime. 

In our study we do not consider finite N corrections to the density of states p 
of the interaction matrix which is taken to be deterministic and does not vary from 
sample to sample. The underlying phase transition in these models is similar to Bose 
Einstein condensation where at low temperatures the system develops a macroscopic 
condensation onto the eigenvector corresponding to the largest eigenvalue p of the 
interaction matrix. In other words the system develops a macroscopic magnetisation in 
the direction of this eigenvector. When the largest eigenvalue of the interaction matrix is 
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bounded these models can present a finite temperature continuous phase transition [TT] - 
[18] . Whether or not this transition is realised depends on how the density of states p 
vanishes at po- In cases with a density of states with tails that extend to infinity there 
is no finite T c , see e.g. [20], [21] . Here we restrict out attention to density of states where 
the maximal eigenvalue is bounded. In the high T phase the system is paramagnetic 
(or liquid) while, as mentioned above, in the low T regime it is ordered via a Bose- 
Einstein condensation mechanism. A finite magnetic field may or may not kill the static 
transition depending on the decay of the density of states close to the edge at po. If the 
density of states behaves as p(p) ~ (p — p) a about the edge at p there is only a finite 
temperature transition if a > 0. In the presence of an external field, this transition is 
killed if a < 1 but there is still a transition if a > 1. In much of our study we will 
concentrate on the regime where a G (0,1). The divergence of all the susceptibilities 
studied here, and other critical exponents, can all be expressed in terms of the exponent 
a. 

From the point of view of the dynamics of these models, a typical initial condition 
does not have a macroscopic overlap with the condensed low temperature equilibrium 
configurations. Consequently the low temperature dynamics is slow and the equilibrium 
condensation is not reached in finite times with respect to N. The relaxation occurs 
out of equilibrium and correlation and linear response functions age with similar scaling 
forms to ferromagnetic domain growth [13] [Hj. In the case where a finite magnetic field 
is applied, and when a < 1 i.e. the field kills the static phase transition, a charactersitic 
time-scale t*(h) is introcuced that separates a transient non-equilibrium regime from 
the final approach to equilibrium in the disordered phase [T5] . 

1.1. Summary of results 

We show that in case (i) the static linear susceptibility does not have a cusp at T c 
and obeys a Curie law down to T = 0. The first static non-linear susceptibility (% 3 ) 
however diverges linearly with the size of the system at all T < T c . In the limit (ii) 
the static linear susceptibility exhibits a cusp and the behaviour of the static non-linear 
susceptibility depends explicitly on the decay of the density of states of the interaction 
matrix at its upper edge. For a density of states decaying as the power law (p — po) a , 
Xs diverges if a < 1/2, vanishes if a > 1/2 and takes a finite value if a = 1/2. 

In a dynamic study, using a Langevin stochastic evolution for the continuous spins, 
we elucidate the approach of the linear and first non-linear susceptibilities to their 
asymptotic static limits in the more interesting low-temperature phase. In short, we find 
that using the first order of limits the linear susceptibility, Xi, approaches its asymptotic 
static value with an exponential decay of characteristic time r = N(/3 — /3 C ). In the 
second case xi achieves its asymptotic value with a power law decay t~ a but the non- 
linear susceptibility X3 diverges or vanishes, depending on a being larger or smaller than 
1/2, with a power law t l ~ 2a - in agreement with the static calculation. In the case (i), 
where the zero field limit is taken first, X3 approaches its asymptotic value exponentially 
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with the same characteristic time as the linear susceptibility, r = N{(5 — (3 C ). 

We stress the fact that the distinction between the two limiting procedures (i) and 
(ii) is especially important in numerical simulations (and experiments). If one takes the 
first order of limits no special feature at T c is observed in the linear susceptibility. One 
needs to use a sufficiently large field and fall into the second case, to see the critical 
behaviour of the linear susceptibility. 

The difference between these two limiting procedures has been discussed in a 
number of papers. Within the droplet model Fisher and Huse derive scaling forms 
for the susceptibility in both cases [22]. More recently, Yoshino and Rizzo |23j studied 
the linear response in mesoscopic disordered models with one step replica symmetry 
breaking using the Thouless-Anderson-Palmer approach and found stepwise signals (as 
anticipated by Kirkpatrick and Young [24] and by Young, Bray, and Moore |25j in their 
studies of the Sherrington-Kirkpatrick spin-glass). They also discussed the importance 
of considering the two limiting procedures (i) and (ii). 

In experiments and numerical simulations, linear and non-linear susceptibilities are 
usually obtained as functions of (higher order) correlations functions by virtue of the 
fluctuation-dissipation relations (in equilibrium) or its extensions (out of equilibrium). 
Some fluctuation-dissipation relations linking non-linear susceptibilities and correlation 
functions in the absence of a field, in glassy systems out of equilibrium, can be found 
in [201 [26]. Numerical recipes to compute the linear susceptibility using fluctuation 
dissipation relations, also valid out of equilibrium, have been proposed by a number of 
authors [271 EEJ EH]- All these expressions correspond to the order of limits (i). 

Having presented our main results above, the body of the paper presents technical 
details of their derivation. In Sect. ([2]) we introduce the model and the notation. In 
Sect. ([3]) we present the free-energy density and the main observable we shall study, the 
magnetisation density. Sect. (j4|) is devoted to the analysis of the Langevin dynamics 
and the derivation of the dynamic susceptibilities. Finally in Sect. (jSJ) we further discuss 
our results and present our conclusions. 

2. The model 

The Hamiltonian considered is for a system N continuous spins 



with the spherical constraint £\ Sf = N. We denote the fixed, or deterministic, density 
of states of the interaction matrix J by 



In what follows we consider the generalised random orthogonal class of models with 




A 



J = O JO 
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where = indicates statistically identical and O is a random rotation chosen with the Haar 
measure. The behaviour of the statics of this class of generalised random orthogonal 
models with Ising spins has been studied in [30] and the first such models were studied 
in [31]. 

We assume that the density of states p(p) has an edge at po and we use density of 
states p that admit the following expansion about p = p 

oo 

p(p)=Y,Cn(V0-Hr +n ■ (1) 

n=0 

For the integrability of p(p) we must have that a > — 1. 

We pay special attention to the The Gaussian ensemble where the elements of the 
(symmetric) Jij are Gaussian of zero mean and variance 1/N. In this case we find the 
Wigner Semi-Circle Law for the density of states 

clearly a = 1/2, p = 2 and Co = 1/tt. This model is statistically invariant with respect 
to such random rotations or transformations. However it does not strictly belong to 
this class of models, as the density of states fluctuates from sample to sample. However 
these sample to sample fluctuations are unimportant for the computation of extensive 
thermodynamics quantities. This model is then equivalent, in the thermodynamic 
limit to the spherical Sherrington-Kirkpatrick disordered model studied in a number 
of publications [TT]-[T9]. 

The spherical ferro magnet in d dimensions, and its continuum version, the O(N) 
model, can also be included in this family of models but we shall not discuss them in 
detail here. From the calculational point of view it is practical to consider an applied 
external field h of the form 

hi = h^ , 

where cij are independent Gaussian random variables of mean zero and variance one. 
This means that the distribution of h is also invariant under random rotations, i.e. 

h = Oh. 

Let us note that h 2 = h 2 N where the over-line indicates averaging over realisations of 
the field h. Consequently the central limit theorem tells us that for large N we have 
h 2 = h 2 N + 0(^/N), this means that after a suitable rotation hi = h + 0(l/y/N), i.e. 
up to 1/yN corrections, there is a basis where h is a uniform vector. The magnetisation 
induced in the direction of the field h is given by 

i 

and for the reasons given above, in the large iV limit m will be a function of h. 
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Pf(z) = \j d ^ p^) ln ( z ~ ^ ~ J d ^ 



3. The statics 

3.1. The free- energy density 

Standard arguments lead to a variational expression for the free-energy (up to constant 
terms irrelevant for our calculations) per spin [TT] . 

f3h 2 f J p(/i) _ f3z 
z — fi 2 

where z is the Lagrange multiplier introduced to impose the spherical constraint. The 
free-energy per site / is obtained via a saddle-point calculation over z as 

/ = min 2 f(z) , 

and the corresponding saddle point equation is 

Here we have used the short hand notation 
<(...)>^ K lim 

to represent the average over eigenvalues. 

3.2. The phase transition 

Whether or not this model exhibits a phase transition depends on the density of states 
of the interaction matrix J^. The constraint equation may be written as 

F(z,P,h)=P (2) 

with 

F{z,P,h) = ((-±- ))+fih* « 7 -L-». (3) 

Z [A \Z \l) 

In this analysis one is required to choose a solution z > /i where /i is the largest 
eigenvalue of the interaction matrix. If F(z, f3, fi) diverges as z — > fiQ then the constraint 
equation can always be satisfied in a continuous manner and there will be no finite 
temperature phase transition. 

In the case h = there is a transition if a > 0. In the case ft ^ there is a 
transition if a > 1. Thus the presence of a finite field for a G (0, 1) kills the transition. 

For models where a G (0, 1) the field has a very singular effect on the 
thermodynamics. In this range of a and in the absence of a field there is a critical 
temperature T c — 1//3 C defined by 

/?c=((- J — ))• (4) 

For (3 < p c , the Lagrange multiplier at zero field z Q varies continuously with (3 and is 
the solution to 

((- 1 -))=/5, (5) 

Z — fi 
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whereas for (3 > (3 C we have z = /x up to 1/N corrections and the spherical constraint 
is satisfied by a macroscopic condensation onto the eigenvector corresponding to the 
eigenvalue fio in a manner analogous to Bose-Einstein condensation. More specifically, 
we decompose ((■■■)) into two parts, the first being over the largest eigenvalue fi and 
the second being over the remaining eigenvalues, which can be treated as a continuum, 
which we denote by ((. . .)) c . In this notation the constraint equation becomes 

11 1 

+ «- -»c = /9- 



N z - fi z - fi 
To leading order in 1/N the solution to this equation is 

3.3. The susceptibilities 

The magnetisation in the direction of the applied field is given by 

m(A) = -I = ,l « iV> • (7) 

In order to fully determine the magnetisation as a function of the applied field strength 
we need to know how z varies as a function of h. From the form of the free-energy we 
see that z must be an even function of h and thus 

oo 

z(h) = J2znh 2n , (8) 

n=0 

at least at high temperatures. The magnetisation then has an expansion in terms of the 
generalised susceptibilities 

oo 

»*(*) = E^V"" < 9 > 

71=0 V ' 

or, equivalently, the n-th order susceptibility is just 

Qn+lf 



d n m(h) 



h=0 un 



(10) 

h=0 



dh r - 

The cases (3 > f3 c must be considered separately as the starting point for the 
expansion, the value of z in equation (jSJ), is different. In both cases we expand 
equation (j5J) in powers of h using the expansion of equation (jSJ) for z. 

3.3.1. High temperatures, T > T c (h = 0). First we consider the high temperature 
regime (3 < (3 C . Using the variables g n defined by 



we find that 



z {h) = z + {3h 2 - P 2 ^ h 4 + 2(3 3 ^ h 6 

92 92 



92 
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and 

m{h) =(3h-Pg 2 h* + 2P% h 5 - ^ + h 7 

92 

+ 2 ^4 (%3^4 + 2^5) ^ 9 + 0(/ill) (13) 
#2 

Note that from the zero field constraint equation for zq one has that go = (3. The first 
three terms look particularly simple, suggesting a simple pattern which however breaks 
down at the next order. From this expression one can read the n-th order susceptibility 
Xni see equations (jHJ) and ffTUj) . Xi — Pi X3 — — 6/3g 2 and so on and so forth. The linear 
susceptibility has a Curie- Weiss behaviour all the way up to T = T c . As one approaches 
T c the Lagrange multiplier z approaches /io and so we write 

z ~ [i + 5z (14) 

where 5z is small with respect to /io- 

The integrals in the definition of the parameters g n , with n > 1 are dominated by 
the divergence of the denominator at the edge /i ~ /io for all a < 1 and therefore 

OO Q /-OO c 

' --a+l-n I 



g n ^c / de — - = c (<^) Q+i n / 



o (e + fc) n io (1 + w)"' 

The scaling of 8z with T — T c is obtained from the analysis of equation (J5]). Introducing 
( ED in © one finds 

« 7 T r % T»=&-0 (15) 

(/i + d2 - /i)(/i - A 4 ) 

which for 5z small gives 

c /(«)(^) Q ^(/? c -/?) (16) 

where 

U a ~ l 71 



1(a) = du = - — r — — . (17) 

' 1 + u sin[7r(l-a)] 



o 



This thus yields 



c I{a) 



and consequently 



c / du 







1 + u 



P~Pc 



-| l+(l-n)/« 

~ (T-T c ) 1+(1 - n)/a . (19) 



From this we find that close to T c , all the non-linear susceptibilities behave as 

This expression diverges as soon as n > l+2a. This means that for a G (0, 1) all Xn with 
n > 3 diverge. Note that the order of limits (i) N —>■ oo h — > or (ii) h — > N — > oo is 
irrelevant at high temperatures. 
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3.3.2. Low temperatures, T < T c (h = 0). At low temperatures we have to consider 
separately the two thermodynamic limits mentioned in the introduction. 

(i) N -f oo h -> 

In this case the solution to the constraint equation for z is given by equation ([6]) 
therefore zq is always strictly greater than /z and the expansion given by equation ( {TBI 
is still valid as the expressions for the g n are still finite and we find that 

g n » iV"" 1 ^ - &)» + (( , \ » c 

^0 + N (J3-p e ) ~ V 

« N n ~\(3 - (3 c ) n + O ({N(i3 - (3 c )) n ~ l - a ) (20) 
As an example, the first non-linear susceptibility is given by 

X 3 = ~W - Pc) 2 N + OiN 1 ' ) . (21) 

(ii) h->0 N -y oo 

When the thermodynamic limit is taken in the presence of an applied field, the 
expansion for /3 > j3 c is carried out with zq = /io but we do not assume an analytic 
expansion for z and thus we write 

z = no + s(h) , 

where s(h) — > as h — * 0. The constraint equation now reads 

((t 77Tt))-^ 2 ((t 1 )) =f 3 c -f3. (22) 

Expanding equation fT22|) to leading order in s(h) and examining the limit /i->0we see 
that to lowest order in h the function s(h) is given by 



s(h) = h 



2 

1—a 



f3c al(a) 



((3 - (3 C ) 

where Co is defined by the expansion of the density of states equation (CD) about fi . The 
key point in this derivation is that the most divergent contribution to the integrals in 
equation ( 1221) as s — >• come from the first term in the expansion of p about /io- This 
most divergent term may be extracted as follows. If //* is the minimal eigenvector of J, 
a typical term to evaluate is 



(fiQ — fi + s) J 1 + u J 1 + u 

as s — > 0. Substituting this result into the expression for the density of states we find 
that to the first two lowest order terms in the expansion 

_ _ a. 

Pcoal(a) 



m(h) = (3 c h — c I(a) 



l + a 



. (23) 



The linear susceptibility is Xi — Pc leading to the appearance of a cusp at T c (where 
X\ passes from the Curie- Weiss law to becoming a constant). 
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l+a 
1-a 



Interestingly the case where X3 exists for T < T c corresponds to the case where 
3 i.e. a = 1/2 which is the Gaussian spherical SK model [H]-[T9]. In the case 
where a < 1/2 we see that X3 is infinite for T < T c , but if a > 1/2 it vanishes, and 
so X3 = 0. The spherical SK model thus turns out to be the marginal case for the 
non-linear susceptibility X3- We find that for T <T C the first non-linear susceptibility 
behaves as 

X3 = ~^cl (-J-) (24) 



while approaching T c from the high temperature region \3 diverges as 

RR T( \( C J ( a )\ o 2 2 ( P 

X3 = -6(JacoI{a) [ - ) = -3n c 



o 2 2 
X3 = -37T C 



which along with equation (|24|) means that in the critical region we may write 

g 

i. e the coefficients of the power law divergence in X3 are the same on each side of the 
transition. 

The Gaussian spherical SK model can be treated non perturbatively and solved 
directly as we show in the next Subsection. 

3. 4- The Gaussian case 

Let us focus here one the case a = 1/2 in which the interaction matrix has Gaussian 
distributed elements and the density states is given by the Wigner Semi-Circle law, 
equation (j2j). This is the spherical SK model. In this case T c = J. The integrals over \i 
in the constraint equation can now be carried out explicitly yielding 



Ph 



from where one easily obtains z(T, h) for all values of the parameters T and h with 
no need to use a perturbative expansion. The numerical representation of the solution 
as a function of h for three values of the temperature, T = 0.5, 1, 1.5, is displayed in 
Fig. [T]-left (with solid lines). These results are compared to the first order terms in the 
series expansion valid for h — > in Fig. [T]-right (with thin lines) . We see that the high 
(T = 1.5), critical (T = 1) and low (T = 0.5) temperature curves are indeed very close 
to the h 2 (high and critical Ts) and h 4 (low T) predictions of Sects. I33TT1 and 13.3.21 when 
h is relatively small. 

The magnetisation in the direction of the field, equation (J7J), reads in this case 

m = — z — \/ z 2 — 4 

The non-linear susceptibilities can also be worked out in detail and agree with the 
results of the previous section in the case where c = 
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Figure 1. Left: the Lagrange multiplier z as a function of the strength of the applied 
field h. Right: check of the power law approach to z(h = 0) for h — > (h 2 at high and 
critical temperatures and /i 2 /( 1_Q ) = h 4 at low temperatures). 



3.5. Phenomenology 

With the aim of testing the relevance of this family of very of simple models to describe 
real spin-glass systems, we compare the critical exponents computed here with the ones 
measured experimentally by Levy and Ogielsky [H [2] , obtained numerically by a number 
of authors [32], and proposed by Fisher and Huse on the basis of the droplet model [22J. 
In this way we try to find an optimal value of a to match the experimental results. 



3.5.1. Comparison with experimental and numerical results. Levy and Ogielsky [HE] 
measured ac non-linear susceptibilities in a dilute AgMn alloy with the characteristics of 
a Heisenberg spin-glass in three dimensions. In their experimental studies they identified 
a finite critical temperature and studied the critical singularities of the nonlinear 
susceptibilities in the static limit. 

Above T c Levy and Ogielsky found 

X3~(^^) with 7 = 2.3 ±0.15, (25) 



and 



— ~ ' -— ~ ' xl +t3h with /3~ 0.9 ±0.2. (26) 



X3 X5 

In the spherical disordered models we found 

-(l-q) 

X 3 ~ (T - T c ) a 



so, comparison with ( 1251) implies 

a ~ 1/3.3 ~ 0.3 . 
In addition we found, 

_^ ^ ^ (T - T c )- l/a ~ x\ /(1 ~ a) 

X3 X5 
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and thus the hyperscaling relation (l2"B"j) gives 

i +/V7 „ x i/(i- a) with ^ ^ = 1 

a 

for all spherical models independently of a for a G (0, 1). Note that /3 = 1 is consistent 
with the analytic behaviour of the order parameter q ea in the low temperature phase. 

It is interesting to note that, as summarised in a recent review article [32], most 
isotropic spin- glasses have ft ~ 0.9 — 1.1 and 7 ~ 1.9 — 2.3. The first value is the one 
we found for all spherical models, the second one implies a ~ 0.3. A scenario including 
a decoupling of spin and chiral order in Heisenberg spin-glasses has been proposed by 
Kawamura [33J. 

As for Ising spin-glasses, both experiments and simulations point to a larger value 
of 7. Kawashima and Rieger |32j stress that it is now well established that there is a 
conventional second order finite temperature transition with a diverging xsg [34"l 135] 
and quote, basically, ft ~ 0.5 and 7 ~ 4 for these 'easy-axis' systems. Daboul, Chang 
and Aharony [36] estimated 7 for the Ising spin-glass model on a hypercubic lattice 
in d > 4 with different distributions of the coupling strengths using high temperature 
expansions. The values they find for higher dimensions also suggest a rather high 7 in 
d= 3. 

We then conclude that, as expected, the spherical model is more adequate to 
describe the high temperature critical behaviour of isotropic rather than Ising-like 
systems. 

Below T c the dynamics are so slow that Levy and Ogielsky could not identify a 
static limit in zero applied field. Aging effects come into play [37] and one has to 
analyse experimental, as well as numerical, data very carefully. Levy's data in a finite 
field are consistent with 7' = 7, with 7' defined from X3 which is finite below T c in the 
experimental case. The value a ~ 0.3 that we extracted from the high T analysis leads, 
however, to a diverging % 3 both in the zero applied field limitbelow T c , see equation (1231) 
obtained with (ii) h — > iV — > 00, and in the opposite order of limits, see equation (12TT) . 
obtained with (i) N —>■ 00 h — > 0. Not surprisingly, spherical disordered models cannot 
capture all details of real spin-glasses. 

3.5.2. T < T c , comparison with the droplet model. We now compare the critical 
behaviour of the spherical disordered models to that of the droplet theory of spin- 
glasses |22j. We shall distinguish the predictions of the latter for Ising and continuous 
spins. 

In the Ising spin-glass phase, T < T c , and in the limit (ii) h — > after the 
thermodynamic limit iV — > 00, Fisher and Huse propose 

m{h) ~h + Ah d ^ d - 2 ^ 

while we have 

m(h) ~h + Ah {l+a)/{1 - a) . 
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An equivalence between the two implies 

a = ^- < 28) 

In the reversed order of limits (i) iV — > oo h — > 0, Fisher and Huse have 

and since we find X3 ~ ^ a comparison leads to 

0(1 + 0) = O # = or = -l. (29) 

Using equation (I2"8"j) the first option, 6* = 0, yields a = 0. Note that 9 = is usually 
associated to the replica symmetry breaking scenario and it has been found numerically 
in d = 2 for J y = ±1 [38]. 

Jonsson e£ a/'s experimental results for the dynamic relaxation of the AgMn 
Heisenberg spin-glass compound analysed with the (slightly modified) droplet scaling 
imply 9 ~ 1 [39]. If we still use the relation (1281) and fix 9 = 1 we then conclude 
a = 1/2 (setting d = 3). Interestingly enough, a = 1/2 corresponds to the spherical 
SK model but also the spherical ferromagnet in the continuum limit (the Laplacian in 
three dimensions leads to a density of states approaching the edge with this power). 

4. Langevin Dynamics 

In this section we compute the temporal behaviour of the magnetisation in the direction 
of the applied field as a function of time for a system quenched from infinite temperature 
at t — 0. The dynamics we study is Langevin dynamics as is the case in most previous 
studies of the dynamics of the spherical SK model [TT]-[T9] 

In the basis where the matrix J is diagonal the stochastic evolution equations 
describing the Langevin dynamics of the system are 
ds 

= (/i- z)s fl + ha^ + r]^. (30) 

In the basis of the eigenvalues the are again uncorrelated and of unit variance. The 
white noise terms have correlation function 

(v,(t) V , / (t')) = 2T5^5(t-t'). 

The term 2 is a dynamical Lagrange multiplier which enforces the spherical constraint 
and which must be calculated self-consistent ly. The solution to equation (|30~!) is 

/ \ / n \ exp(ut) exp(ut) f . 
s»(t) = s M (0) — + ds (hcr^ + rj^s)) exp(-//s) T{s) 

where T(s) = exp ^ Jq ds z(s)j. Now assuming that the initial conditions are such that 
they are uncorrelated with the applied field and also assuming that (s M (0)s At '(0)) = 5^> 
we obtain the following equation for the magnetisation in the direction of the applied 
field: 

1 h f l 

m(t) = j^^^ns^t) = ^-y ds ((exp[/i(t- s)]))T(s) 
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(the field is applied at the preparation time t = and subsequently kept fixed). The 
self-consistent equation for V is 

T 2 (t) = (( exp(2/ii) )) + 2T f ds (( exp [2/i(t - s)] )) T 2 (s) 

Jo 

+ h 2 [ dsds' (( exp \n(2t - s - s')} )) r(s)r(s') . (31) 
Jo 

We restrict out attention to the dynamical behaviour of just the linear and first non- 
linear susceptibilities. The above equations are thus solved perturbatively to 0(fi 3 ) to 
give 

m(t) = h f ((exp [^ -*)]>> ^ 
Jo L o\t) 

+ h ^p s{Mt - sm ^)^pEM , (32) 

where T obeys 

rjj(f) = (( exp(2pt) )) + 2T [ ds ((exp [2/i(t - s)} )) T 2 (s) , 

Jo 

and is given by 

2r„(*)r 1 (*) = / dsds' ((exp [fi(2t - s - s')\ » r (s)r (s') 

Jo 

+ 4T f ds {{exv[2yL{t- s))))T Q {s)T x {s) . 
Jo 



(33) 



(34) 



The dynamical linear susceptibility is then 



Xi(t)= fdsiiexpi^t-s)]))^, (35) 
Jo L o{ t ) 

and we may write the dynamical nonlinear susceptibility as 

xsW = -6^|| + f^y / rfs ((exp ^ - s)] )} Tl{s) • (36) 

4-1. Low temperatures 

In [H] the solution of equation (1331 for a general p(/i) in the ageing regime T <T C was 
found, this result can be compactly written as 

rg(t) = y d f j,q(n)exp(2 fJ ,t), (37) 

with g(/i) given by 

'<"> = ^f&F (38) 

and 

= pJdX^-, (39) 
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where P denotes the principal part. For the sake of completeness we re-derive the 
result equation fl38|) in a new more direct way. First if we assume the representation 
equation ( 1371) we find 



dp, p(p) exp(2/ii) = J dp q(p) exp(2/it) 

, rp [ i A , / a / \ [ exp(2/i / t) -exp(2/it) ] 

+ T j dpdp p(p ) q{p) _ . (40) 

Notice that the apparent singularity in the second integral on the right hand side at 
p = pi is not really present and we can replace the integral by its principal part. 
Equating the coefficients of exp(2pt) in the above equation now yields 

q(p) = p(p) +T X (p) q(p) + T p{p)P [ dp' . (41) 



p- p' 



The Laplace transform of equation (1331) reads 

rg(p)= p ~ 2At 



2T/d//- 



p-2/i 

where 



rg(p) = f°dt exp(-p*)rg(t) = /rf/i . (42) 

The above now implies that 

Pi«m. ^J^_ = _XM . (43 ) 
' /i-/ i-TPjdp'^ i-^x(a0 

Using this result for the last term in equation (j4"Ii) we obtain equation ( 1381 ). 



The Gamma function 

Let us analyse the asymptotic behaviour of T in the two limits. At late times the 
dominant contribution to rj^t) comes from around p = po. Expanding about this point 
we find 

where we have used equation (CQ) for the density of states p at the edge of the spectrum. 

In what follows without loss of generality we restrict ourselves to the case where 
Po = which can be achieved by a constant shift in the energy by using the interaction 
matrix J' = J — p I. From the definition of xif 1 ) m equation (1551) and equation (j3J we 
find 

2 c T(l + a) 

Lo[t) (2i)i+«(l- (44j 

where T in the above is the standard gamma function defined as 

POO 

T(x) = dt r- 1 exp(-t). 
Jo 
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We now compute the large time behaviour of V in the region T < T c . Notice that 
the Laplace transform of T at small p is dominated by the large t behaviour of r (t), 
for p small and a < 1 we have that 



r (] (/')~.-! / iltrx\A-pt)t~^ ~ A V ( — ^ ) 



(45) 



where A = ^/c T(l + a)2-( 1+Q )/ 2 (l - T/T c ). 

In the case in which we keep the N dependence, useful to study case (i), we find 

1 1 



r (p) 



y/N{l-T/T c ) V-W(^r]T c 



(46) 



^.S 1 . The linear susceptibility 

Equation (I35p can now be written as 

K{t) 



with 



Xiit) 



K(t) 



r (*) ' 



als ((exp [//(t - s)])) r Q (s) . 



The Laplace transform of is given by 

*(p) = ((A-)) rJ(p) • 



We may thus write 
K{p) 



p-fi 



((--))+ P ((^—, 
p p(p - p) 



r (p) 



T c p{p - p) . 



r (p) • 



fii) /i -> iV -> oo 



(47) 



The term ((l/(p(p — //)))) diverges as p — > for a < 1; the small p behaviour is 
thus dominated by the region around p = 0. This gives for small p 



«- 



;)) ~ -c / <fe 



(48) 



>(p-/i)' 

where 1(a) is as defined by equation ( flTj) . From the above and equation ( 1451) we thus 
find that for small p 

tf(p) « - co/(a)Ar(^-)p— . 

Asymptotically inverting the Laplace transform we obtain 



r (t) c /(a)Ar( 
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which finally yields for large t 



Xi(t) 



l-3a-> 



(49) 



We thus see that Xi(t) decays to its low temperature equilibrium value with a power law 
t~ a . Interestingly the coefficient of this decaying term is negative for a < 1/3, meaning 
that Xi(t) achieves its equilibrium value from below, whereas for a > 1/3 the coefficient 
is positive and thus Xi{t) achieves its equilibrium value from above. 

(i) iV -> oo h -»• 

In this limit one can analyse K(p) in ( 14"T1) keeping the 1/N contributions. One finds 

I-fl-I 



Xi(t) 



g N(f3-/3 C ) 



Note that if we set t/N <C 1 we recover \i ~ A- 
4-4- The non-linear susceptibility 

We now turn to the results for the nonlinear susceptibility x 3 . If we define Q(t) = 
T (t)T 1 (t), from equation (I34|) we find that the Laplace transform of Q obeys 

where 

1 f* 

/(*) = o / dsds ' ^ exp t 2/it ~ ^ s ~ V s 1 )) F o( s ) r o(s') ■ 
(ii) /i -> iV -> oo 

Let us now focus on this order of limits. Making the substitution w = —/it and 
s = vt we find 



m= \L 



1 r"** f-w\ dw 

p 





t 



1 -.2 

tdu exp [— w(l — v)] T (vt) 



We now use the asymptotic form of r in equation ( 1441) to find, for large t, 



f(t) 



A 2 c 
2t 2a 



dw w a 



dv 



exp [— w(l — v) 

l + ct 
V 2 



Thus f(t) = B'/t 2a for large t and one may verify that the constant B' is finite for 
a < 1. Consequently, we obtain 

r x (t) = st^ 1 - 30 ) . 

The small p behaviour of the Laplace transform of Ti is thus given by 

3 - ^ -*(l-a) 



r^^rl-fi-ajjp 
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We now rearrange the result in equation (I55T) as 



Xs(t) = -6 



ri(t) 
r (t) L 



6 



L(t)_ 

r (t) 



where 



The Laplace transform of L is given by 

1 ,,1 



L(p) 



» r 1 ( P ) = - P r 1 (p)« 



(50) 



For small p using the asymptotic result for Ti and equation fHHl) we find 



Up) 



c BI(a)T(-(l 



a) p 



which implies that the late time behaviour of L is 

c B 1(a) T (f(l - a)) 



L(t) 



r(l 



fa) t2 a -2 



We can compute the other contribution to X3 using the asymptotic result for xi 
equation (09]) to obtain the result 



6c B 1(a) 
A 



r (|(i -a)) r( 



1-a' 



r(l 



7« 



•l-3a^ 



l-2a 



2"V V 2 

We therefore see that when a < 1/2 the nonlinear susceptibility diverges as x^) ~ 
Ct 1 - 2 *. When a > 1/2 xs(t) decays to zero as l/t 2a 1 . These dynamical results are of 
course in agreement with the static calculations carried out earlier. The coefficient of 
this term is determined by the sign of that in the square brackets, the first prefactor 
being negative. Numerical evaluation of the factor in square brackets confirms that it 
is positive for a G (0, 1) and thus the dynamical nonlinear susceptibility is negative. 
Indeed for T > T c the static value is negative and the dynamical calculation confirms 
the divergence to an infinite negative value for a < 1/2. 



fi) N -> oo h 







In order to compute X3(t) we need to compute T (t) and Ti(t) keeping the correction 
to the leading exponential in time terms. From equations ()34p and (1461) we find 

1 



r (t) 



Ti(t) 



N(l- T/T c 



e N(fS-0 c ) 



1 



2T 2 



y/N(l-T/T c ) {1 + 2t/[N(P - p e )]} e^^y 



(51) 



(52) 



Using equation ( 1361) it can be verified that the asymptotic limit of Xs(t) * s the static 
value ()2ip up to a correction term that decays exponentially as e~ N (P-0c) . 
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We have studied in detail the linear and first non-linear susceptibilities of generalised 
random orthogonal model spherical spin glasses. Their physics is completely determined 
by the density of states of the two body interaction matrix. In particular the exponent, 
a describing how the density of states vanishes at the upper edge, determines completely 
the critical behaviour at the phase transition and the dynamical evolution of Xi{t) an d 
Xsit) in the limit (ii) considered here where the limit h — > is taken after the limit 
N — > oo in the computations. With respect to the Gaussian p = 2 spin glass model 
we see that the existence of the parameter a gives us the possibility to carry out a 
more meaningful comparison between the model and experimental and droplet scaling 
theories. An interesting aspect of this work is that it clearly demonstrates the possibility 
that X3 m &y appear to be finite for certain classes of model (with a > 1/2) if the applied 
fields used to carry out the measurements place us in the regime of limits (ii). However 
it is in this same region where the cusp in the linear susceptibility exists. These results, 
though for a somewhat idealised mean field model, could well have some bearing on 
the interpretation of susceptibility and magnetisation measurements in disordered and 
frustrated spin systems [10]. Finally we have mentioned that the models considered here 
will still exhibit a finite temperature transition in the presence of an external magnetic 
field if a > 1. In this case it is the higher oder (n > 3) susceptibilities which will diverge 
and it will be interesting to study, in particular, the dynamics of these models. 
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